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Characterization of the interaction between dissociated air species and 
thermal protection surfaces is prerequisite for theoretical prediction of surface 
catalycity phenomena under re-entry conditions. Frozen orbital methodol- 
ogy and software necessary for the ab-initio quantum mechanical calculation 
of accurate gas-surface interaction potentials has been developed at LaRC 
under the current contract. The codes compute energies, properties and 
wavefunctions for large, spatially extended systems using a bootstrap frozen 
orbital approach. 

The codes were applied to several small linear crystalline test systems. 
Energy shifts caused by variation in the structure of central defect sites are 
accurately predicted when the bulk-local site interaction was small. In other 
instances (more representative of those encountered around recombination 
sites) observed discrepancies with all electron results are due to: 

1. Basis set superposition errors (the bulk basis set is inaccessable to the 
variational treatment and therefore the model results are not subject 
to spurious lowering of the energy due to superposition of the free and 
frozen bases) 

2. Omission of bonding effects at the boundary in the absence of basis 
boundary transformations. These errors can be alleviated by increasing 
the size of the local site and/or by using boundary transformations to 
efficiently repartition the basis. 

3. First order bulk polarization errors. When the defect field is very 
strong, bulk orbitals can be correctly polarized by self-consistently op- 
timizing them in the defect field. 

4. Second order dispersion and charge transfer errors. These errors can 
be minimized by increasing the size of the local site. 
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Boundary bonding and bulk polarization effects had the most serious impact 
on accuracy of the model gas-surface potentials. Further methodological 
work will be directed towards an improved treatment of boundary effects 
and towards minimizing the bulk polarization error by means of parallel 
heterogenous solid calculations. 

An application and thorough test of the method involves prediction of the 
effects of hydration on the electronic structure of the nitrate ion by freezing 
the description of the second solvation shell. Initial studies are described in 
the attached paper “Ab Initio Studies of Nitrate Ion Hydration Complexes". 
The self-consistent field method with accurate basis sets has been applied 
to compute completely optimized equilibrium geometries, vibrational fre- 
quencies, thermochemical parameters, and stable site labilities of complexes 
involving 1, 2, and 3 waters. This study will be submitted to the Journal of 
Chemical Physics. 



Ab Initio Studies of Nitrate Ion Hydration 

Complexes 


Frederick Senese 1 

Old Domtnton University Research Foundation. P. O. Box 6369. 800 W. 

46th St., Norfolk VA 23508 
Donald H. Phillips 

Mail Stop 231, NAS A -Langley Research Center, Hampton, VA 23665 

Abstract 

Ab initio quantum mechanical calculations on small hydration com- 
plexes involving the nitrate anion are reported. The self-consistent 
field method with accurate basis sets has been applied to compute 
completely optimized equilibrium geometries, vibrational frequencies, 
thermochemical parameters, and stable site labilities of complexes in- 
volving 1, 2, and 3 waters. The most stable geometries in the first 
hydration shell involve in- plane waters bridging pairs of nitrate oxy- 
gens with two equal and bent hydrogen bonds. A second extremely 
labile local minimum involves out-of- plane waters with a single hydro- 
gen bond and lies about 2 kcal/mol higher. The potential in the region 
of the second minimum is extremely flat and qualitatively sensitive to 
changes in the basis set; it does not correspond to a true equilibrium 
structure. 


1 Introduction 

Nitrate ion hydration complexes play an important role in elementary mech- 
anisms of cloud nucleation [1, 2], polar stratospheric cloud chemistry [3], 
D-region negative ion chemistry [5], and tropospheric NO r and acid depo- 
sition chemistry. The structure and energetics of nitrate hydration clusters 
is therefore of environmental and aeronomic interest. Thermodynamic prop- 
erties of hydrated nitrate anions have been investigated using high-pressure 

'Present Address: Dept, of Chemistry, Frostburg State University, Frostburg MD 
21532 
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mass spectroscopy by Lee [6, 7], Payzant [8], and Fehsenfeld [5] and cowork- 

ers. • 

Lee and coworkers [7] performed CNDO/2 calculations on one, two, and 

three water nitrate anion hydrates. Their calculations showed bent hy- 
drogen bonds with lengths of 1.38, 1.40, and 1.42 A, respectively, with 
_ o — N = 115°. Configurations with water lying out of the plane of 
the anion were slightly favored over planar complexes. The first ab tnxtto 
studies of nitrate anion hydration complexes were performed by Howell and 
coworkers [12]. Their study treated one-, two, and three-water complexes 
at the 6-31G/SCF level. Geometry optimizations were severely constrained, 
with only planar complexes with water fixed at its experimental equilibrium 
geometry being considered. Single hydrogen bonds were further constrained 
to be linear. The two and three-water complexes were restricted to a few 
high symmetries. No vibrational analyses were performed to verify that 
the geometries were true local minima. A recent SCF study by Shen, Xie. 
Schaefer, and Deakyne[13] reported completely optimized stationary points 
of the monohydrate potential surface using a double-zeta plus polarization 
basis set. Anion diffuse functions were not included in the basis since these 
functions have relatively little effect on computed binding energies; however 
it is generally accepted that diffuse functions are necessary to provide an 
adequate description of anionic charge densities [14] and hydrogen bonded 
complexes. A structure with two hydrogen bonds (with a binding energy 
of 14.9 kcal/mole) and a structure with a single hydrogen bond (with a bind- 
ing energy of 12.3 kcal/mole) were reported. Both structures were planar. 
In a subsequent study of isomers of the conventional nitrate monohydrate 
complex [15], single-point calculations including anion diffuse functions were 
performed at geometries which were optimized without the benefit of this 
refinement. 

In the present study, the dependence of the qualitative and quantitative 
aspects of the SCF potential energy surfaces of singly, doubly, and triply 
hydrated nitrate anions on the quality of the basis set has been investigated. 
Additionally, thermodynamic properties of NO^H^O^ complexes have been 
computed at the SCF level using a valence triple zeta quality basis set with 
both polarization and anion-diffuse functions. Geometry optimizations were 
completely unrestrained. Vibrational analyses using computed first and sec- 
ond energy derivatives were performed for each optimized geometry. Equilib- 
rium concentration ratios of NO^O)* complexes under ionospheric and 
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polar stratospheric conditions are predicted using the ab tnitto thermochem- 
ical parameters with a standard atmosphere model. 

2 Method 

Ab initio self-consistent field calculations and vibrational analyses were per- 
formed on NO3 (HaO) n complexes (with n = 1,2, and 3) using the GAMESS 
program [17]. Three basis sets were used: 

• A Dunning/Hay double zeta plus polarization (DZP) basis set with 
scale factors of 1.2 and 1.15 for hydrogen [16, 17] was used for prelimi- 
nary exploration of the potential energy surfaces. This N, O (10s5pld)/[3s2pld], 
H (4slp)/[2slp] basis set contains 85, 110, and 135 contracted Gaus- 
sian functions for NO^FUO),, calculations with n = 1,2, and 3, re- 
spectively. 

• An extended Pople 6-311G** basis set[18] augmented with anion diffuse 
sp shells (( = 0.0845)[19] centered on each nitrate oxygen was used to 
obtain fully optimized one- and two-water complex geometries. With 
this triple zeta plus polarization plus diffuse basis set (TZP +) there are 
119, 150, and 181 contracted Gaussian functions for the singly, doubly, 
and triply hydrated complexes, respectively. 

• An N, O ( 1 3s8p2d ) / [7s4p2d] , H (8s2p)/[5s2p] van Duijneveldt basis set 
[20, 21] augmented with anion diffuse sp shells (( = 0.06) on the nitrate 
oxygens was used to compute N0 3 (HoO) complex binding energies and 
charge density maps at the TZP+ optimal geometries. This quadruple 
zeta plus double polarization plus anion diffuse (QZ2P+) basis set con- 
tains 189 contracted Gaussian functions for the monohydrate complex. 

The accuracy of SCF level hydrogen bonding energies is limited by the 
ability of the basis set to describe the properties of the individual fragments. 

Errors in fragment dipole moments, static polarizabilities and ionization po- 
tentials will propagate to errors in the electrostatic, polarization, and charge- 
transfer contributions to the hydrogen bond energies, respectively. [22] Ener- 
gies and structures computed for the isolated nitrate and water fragments are 
compared with those obtained by experiment and previous SCF calculations 
in Table 1. 
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Table 1: Optimized geometries and properties oj isolated H?0 and NOj frag- 
ments. Distances are m A, angles are in degrees, energies are in a. u., and 
dipole moments are m Debyes. 


Fragment Source 

Optimized Geometry 

Total Energy 

Dipole Moment 

H,0 

-Roh ZHOH 



1 

0.940 106.4° 

-76.06370 

2.00 

2 

0.941 106.9° 

-76.060411 

2.25 

3 

0.941 104.4° 

-76.04709 

2.14 

4 

0.941 107.1° 

-76.05604 

2.25 

6 

0.944 106.7° 

-76.04683 

2.19 

7 

0.95781 104.4776 


1.8546 

NO3 

R NO 



1 

1.222 

-279.04831 

0.0 

O 

Ad 

1.224 

-279.03967 

0.0 

3 

1.222 

-279.01063 

0.0 

4 

1.223 

-279.02533 

0.0 

5 

1.26 

-278.7676 

0.0 

6 

1.231 

-278.99006 

0.0 


1. This work; SCF with QZ2P+ basis. 


2. This work; SCF with a basis that is the same as the QZ2P+ basis but with 
only one d polarization function on 0, N and only one p polarization function 

on H. 

3. This work; SCF with TZP+ basis. The SCF/TZP+ vertical ionization po- 
tentials for water and nitrate anion are 0.4011 and 0.1645 hartrees, respec- 
tively. 

4. This work; SCF with a Dunning triple split valence [23] plus polarization 
plus anion diffuse sp shells on nitrate oxygens. 

5. See Reference [12], SCF with Pople 6-31G basis set. 

6. This work; SCF with DZP basis. The SCF/DZP ionization potentials for 
water and nitrate anion are 0.4023 and 0.1590, respectively. 

7. Experimental values: see References [25] and [26] for water. 



5 


Figure 1: Optimized geometries , total energies and dipole moments of mono- 
hydrate complexes of NOj . 


©© 




Structure 

E total (a. uj ( DZP1 -355.061144 
(TZP + ) -355.080509 

(QZ2P+) -355.133512 

p D (Debye 1 (DZP) 0.140 
(TZP+) 0.452 
(QZ2P+) 0.368 



lb 

-355.056901 

-355.077085 

1.927 

1.801 


The TZP+ energy and dipole moment of water is only a slight improve- 
ment over that of the DZP basis but is significantly better than that provided 
by the 6-31G basis. The SCF/QZ2P+ basis yields a total energy for water 
that estimated to be 0.0038 a. u. above the SCF limit [24] and overestimates 
the gas phase dipole moment of water by only 8%. 


3 Results and Discussion 

Optimized geometries of the mono, di, and trihydrate complexes are shown 
in Figures 1, 2, and 3, with a comparison of geometric parameters shown in 
Table 2. 

The most stable geometries in the first hydration shell involve in-plane 
waters bridging pairs of nitrate oxygens with two equal and bent hydrogen 



Figure 2: Optimized geometries, total energies and dipole moments of dihy- 
drate complexes of N0 Z ■ 




Structure 

E total (a. u.) (DZP) 
(TZP+) 

p D (Debye) (DZP) 
(TZP+) 


Ila 

-431.129326 

-431.148512 

0.438 

0.210 


Structure Hb 

E to tai(a. u.) (DZP) 

(TZP+) -431.145296 

Pd (Debye) (DZP) 

(TZP+) 1-428 



Structure 
Etotal(a- u.) 

p D (Debye) 


~ lie 

(DZP) 

(TZP+) -431.142154 
(DZP) 

(TZP+) 1 ID 4 



I 


Figure 3: Optimized geometry of the most stable tnhydrate of N0 Z . The 
D 2h geometry Ilia is a local minimum in the DZP basis and a double saddle 
point in the TZP+ basis. The arrows on the m-plane view show nuclear 
displacement vectors corresponding to one of the TZP+ doubly degenerate 
imaginary frequencies. Illb is the local minimum in the TZP-h basis. 



Structure Ilia 

EWa. u.) (DZP) -507.195071 

(TZP + ) -507.214732 
HD (Debye) (DZP) 0.00 

(TZP+) 0.00 




(TZP+) -507.214736 
(TZP+) 0.03 
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bonds. A second energy minimum involves out-of-plane waters with a single 
strong hydrogen bond and lies about 2 kcal/mol higher. Computed total en- 
ergies and thermochemical properties of the hydration complexes (and frag- 
ments) are given in Table 4. The thermochemical properties were computed 
using the ideal gas, rigid rotor, and harmonic normal mode approximations 
at standard temperature and pressure. The larger dipole moment predicted 
by the DZP basis in conjunction with greater basis set superposition error 
leads to slightly higher binding energy estimates than those obtained with 
the more extensive bases. 

The NOj ••• HoO Potential 

The local minimum at la closely resembles the doubly hydrogen 
bonded monohydrate structures reported by Howell et. al. and Shen et. 
al. The pair of symmetric bent hydrogen bonds at the SCF/TZP+ level are 
0.02 A longer and about 0.5° more bent than those computed by the latter 
authors with a double zeta plus polarization basis set. 

Figure 4 shows slices through the potential well at la at a) fixed nitrate 
oxygen to water oxygen distances Roo an d b) at fixed water roll angle 8 
about its dipole. The shape of the roll angle curve can be qualitatively 
explained using a simple point charge model; the planar configuration is 
electrostatically favored. 

Figure 6 shows the effect of complex formation on the QZ2P+ charge 
density distribution for structure la. Polarization of the water is the primary 
effect of complexation on the charge density. Figure 5 shows the effect of 
complex formation on the orbitals of the nitrate and water fragments at the 
QZ2P+ level; the electron donor orbitals are lowered and the proton donor 
orbitals are raised [22]. 

A second TZP+ local minimum (lb) involves out-of-plane waters with a 
single strong hydrogen bond and lies about 2 kcal/mol higher than structure 
la. The potential in the region of lb is extremely flat and quite sensitive 
to basis set changes; for example, in the DZP basis lb collapses to a planar 
structure. Two saddle points were located in this region on the TZP + surface: 

1. Ic has C, symmetry lies at -355.076978 and resembles equilibrium ge- 
ometries reported by Lee et. al. and Shen et. al.. at the CNDO 
and double zeta plus polarization level, respectively. However, at the 
SCF/TZP+ level this structure is a saddle point. The single imaginary 



9 


Table 2: NOJ • • HiO complex stationary point geometries. 

A, angles are in degrees. 

Structure B asis Ru jfo Ri* ^ 2g ~^ 37 

la DZP 1.235 1.235 1.222 2.185 2.185 

TZP+ 1.226 1.226 1.213 2.209 2.209 

lb DZP 1-227 1.240 1.224 4.227 2.006 

TZP+ 1.217 1.231 1.215 4.383 1.987 

Ha DZP 1.239 1.226 1.226 2.209 2.238 

TZP+ 1.230 1.217 1.217 2.216 2.253 

Hb TZP+ 1.221 1.219 1.222 2.241 2.225 

He TZP+ 1.211 1.224 1.226 4.840 2.012 

Ilia DZP 1230 1.230 1.230 2.258 2.258 

IHb TZP+ 1-221 1.221 1-221 2.253 2.265 


Distances are in 


Rb6 

0.949 

0.946 

0.943 

0.941 

0.949 

0.945 

0.945 

0.941 

0.948 

0.945 


R 57 

0.949 _ 

0.946 

0.956 

0.951 

0.948 

0.945 

0.945 

0.950 

0.948 

0.944 


Struct ure Basis ^213 Z214 

Hi DZP 119.9 120.2 _ 

TZP+ 1195 120.3 

lb DZP 119.4 120.8 

TZP+ 119.7 119.6 

Ha DZP 119.8 119.8 

TZP+ H9.7 119.7 
Hb TZP+ 120.2 119.9 

He TZP+ 120.3 120.3 

Ilia DZP 120.0 120.0 

IHb TZP+ 120.0 120.0 


Z126 '137 

111.2 111-2 
111.5 111.5 
115.2 
130.9 
110.9 111.9 
111.0 112.0 
111.1 111.6 
134.8 


'265 

139.2” 

138.9 


139.2 99.6 

138.9 99.7 


149.3 

147.6 

140.1 100.4 

140.4 100.2 

138.2 100.2 

102.2 


103.0 

102.0 


111.5 111.5 138.4 137.6 1006 


Structure 


lb 


Basis 


DZP 

TZP+ 

TZP+ 


0.0 

144.8 

-0.7 


180.0 

166.5 


TZP+ -151.3 -12.2 


'2149 

'1498 

Z.498 10 

13.5 

11.0 

143.9 

151.0 

-12.4 

-10.6 


'9 8 10 


Structure B asis 1 H49 #98 to t 1 = 

£ TZP+ 2.019 0.950 0.941 130.2 146^6 102.2 

lie TZP+ 2.015 0.950 0-941 132.3 146.9 102.3 
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frequency leads to equivalent structures lb above and below the plane 
of the nitrate ion through torsion around the hydrogen bond and the 
adjacent nitrate NO and water OH bonds. Ic lies only 0.07 kcal above 
lb so the motion of water between the equivalent lb structures through 
the anion plane is barrier free. 

2. Structure Id is a saddle point along the minimum-energy path between 
structures la and lb and lies at -355.077059 hartrees, conversion of lb 
to la is therefore also barrier free. 

Figure 7 shows a slice through the potential energy surface corresponding 
to torsion about the hydrogen bond to show the relationship between the 
minima at la and lb and the saddlepoint Id. The flatness of the potential 
in the region around lb precludes its consideration as an equilibrium struc- 
ture”. Note that the minimum lb has a D. of 12.5, 12.1, and 11.3 kcal/mol 
at the SCF/DZP, SCF/TZP+, and SCF/QZ2P+ levels, respectively, Shen’s 
singly hydrogen bonded structure has a De of 12.3 kcal/mol [13]. 

Several higher order stationary points were located. A planar C*2v struc- 
ture with two “bifurcated" hydrogen bonds involving a single nitrate oxygen 
lies at -355.076019 hartrees. A second C iv structure with water lying in the 
anion’s o v plane and bisecting the O-N-O angle lies at -355.075744 hartrees. 
Both structures are double saddle points. 

Table 3 lists vibrational frequencies and infrared intensities computed at 
the TZP+/SCF level for complexes la, Ha, and Illb. The solvated nitrate 
ion in dilute alkali metal nitrate solutions gives rise to observed frequencies at 
1404 cm -1 , 1348 cm -1 , 1049 cm -1 , 825 cm -1 , and 719 cm' 1 [32]. Note that 
energy first and second derivatives are consistently overestimated at this level 
of theory [2 7]; a scaling factor of 0.89 is often applied to calculated frequencies 
in comparisons with experimental spectra. While the anharmonicity of the 
potential surface and inherent deficiencies at the SCF level of theory pre- 
clude quantitative prediction of spectroscopic observables, the calculations 
qualitatively predict many features of aqueous alkali metal nitrate Raman 
and infrared spectra. 

The experimentally observed splitting of the nitrate band in solu- 

tion has been attributed to lowering of the anion symmetry by the ion-water 
interaction [32]; our calculations support this interpretation. The experimen- 
tal splittingfor dilute aqueous alkali metal nitrates is 56 cm -1 : the calculated 
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Figure 5: QZ2P+ molecular orbital energies for complex la. 


A' 2 -.2218 . . 
E" -.2420 -. 

E' -.2975 ••• 


B\ -.2465 

B 2 - 2646 
Aj -.2679 


3 : 

B 2 


-.3196 

-.3207 

-.3359. 


Ai -.4300. 



-.5209' ••• 

-.5345* 


A[ -.6321- ... 


no 3 - 



-.5293 

-.5449 

-.5575. 

-.5580 


f?2 -.5109 


rli -.5963 


,4i - 6577 


no 3 ~ • • • h 2 o 


Bi -.7143 


h 2 o 




1 1 I 1 I I I 

Figure 6: QZ2P+ charge density difference plot for la obtained by subtract- 
ing the isolated fragment one-electron densities from those of the complex. 
Hie solid lines are contours of increased charge density due to complex for- 
mation, while the broken lines repiesent lowered density contours. 
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binding energy (kcal/mol) 
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Figure 7: Monohydrate hydrogen bond torsion potential. Geometries with 
the hydrogen bond torsion angle fixed at the indicated angle with all other 
geometric parameters optimized are marked by diamonds. The solid diamond 
indicates the position of the saddle point Id. 



splitting is 39 cm -1 for the monohydrate la and 38 cm -1 for the dihvdrate 
Ila. However splitting of the u^E') band due to solvent water is predicted 
to be very small (being less than one wavenumber for la and less than 5 
wavenumbers for Ila. This is in accord with Irish’s observation that i/t(E ) 
splittings occur only in concentrated solution and appear to indicate prox- 
imity of the cation [33]. 

The lowest 6. 13, and 18 frequencies in the spectra of complexes la. 
Ila, and Illb, respectively, are associated with rather facile motions about 
hydrogen bond. 


Multiple Hydration Complexes 
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Table 3: Computed vibrational spectra for complexes la. Ha. IHb and in- 
finitely separated water and nitrate fragments. Frequencies are m reciprocal 
centimeters; infrared intensities are in Debye* /amu- A * . 


la 


Ila 


Ilia 

NO 3 

+ HoO 

Freq. 

Int. 

Freq. 

Int. 

Freq. 

Int. 

Freq. 

Int. 

804 

0.0 

805 

0.0 

811 

0.03 

799 

0.0 

804 

0.0 

810 

0.0 

811 

0.03 

799 

0.0 

977 

1.0 

974 

1.1 

973 

1.2 

980 

0.8 

1234 

0.1 

1234 

0.1 

1239 

0.0 

1234 

0.0 

1552 

17.8 

1560 

18.1 

1578 

23.0 

1570 

21.0 

1596 

25.6 

1598 

26.7 

1578 

23.0 

1570 

21.0 

1839 

3.2 

1829 

4.2 

1823 

4.8 

1750 

1.9 



1840 

2.1 

1823 

4.9 







1837 

0.0 



4109 

3.2 

4102 

3.0 

4124 

4.0 

4141 

0.4 



4145 

2.8 

4124 

4.0 







4126 

0.01 



4139 

1.8 

4157 

1.7 

4165 

0.006 

4236 

1.3 



4182 

1.8 

4168 

2.5 







4168 

2.5 




N0 3 - 1 /«(£') 
NOj t'i(E') 
N0 3 - u 2 (A'f) 
NOj ui(A[) 
N0 3 - i*(F) 
NO3 1 *(£') 
H2O bend 


HjO sym. stretch 


HjO asym. stretch 


1 . Hydrogen bond vibrational frequencies (and intensities) for structure la are 
32 cm " 1 (0.0), 87 cm " 1 (0.5), 165 cm " 1 (0.2), 227 cm " 1 (2.3), 235 cm ' 1 
(0.0), and 648 cm -1 (5.9). 

2 . Hydrogen bond vibrational frequencies (and intensities) for structure Ila 
are 22 cm -1 (0.02), 34 cm -1 (0.0), 45 cm -1 (0.3), 108 cm -1 (0.4), 150 cm -1 
(0.01), 165 cm -1 (0.2), 185 cm -1 (0.01), 193 cm -1 (0.001), 213 cm -1 (2.4), 
266 cm -1 (2.6), 604 cm -1 (0.04), and 613 cm -1 (11.6). 

3 . Hydrogen bond vibrational frequencies (and intensities) for structure Ilia 
are: 9 cm -1 (0.01), 19 cm -1 (0.04), 29 cm -1 (0.02), 46 cm -1 (0.3), 47 cm -1 
(0.3), 83 cm -1 (0.0001), 127 cm -1 (0.002), 135 cm -1 (0.0003), 144 cm -1 
(0.2), 150 cm ’ 1 (0.002), 156 cm -1 (0.05), 158 cm -1 (0.002), 189 cm -1 (3.9), 
193 cm ' 1 (3.9), 217 cm -1 (0.003), 569 cm -1 (0.02), 569 cm -1 (0.04), 537 
cm -1 (17.3) 
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In structure Ha the water- water interaction is about +1.2 kcal/mole at 
the TZP+/SCF level (assuming that fragment interactions are pairwise addi- 
tive). The net repulsive interaction leads to slightly weaker hydrogen bonds 
but both waters remain in the plane and the angle between water oxygens 
(with vertex at the nitrate nitrogen) does not open significantly (119.7°). 
Induction associated with geminal hydrogen bonding drives Mulliken and 
Lowdin charges on the doubly hydrogen bound nitrate oxygen 12% and 15% 
higher than those on the singly hydrogen-bound nitrate oxygens, but geminal 
hydrogen bonds are still 0.03 — 0.04 A longer than nongeminal bonds. 

Stationary points at lib and lie involve one and two waters, respectively, 
in highly labile singly hydrogen bonded sites. Again, the potential around 
these sites is extremely flat and lib and lie are not well-defined equilib- 
rium structures despite the fact that they are mathematical minima. At the 
SCF/TZP+ level, D e is 25.4 and 23.4 kcal/mol for lib and Ha, respectively. 
The repulsive interaction between singly and doubly waters in structure lib 
hydrogen bound sites is about 1.0 keal/mol, while repulsion between pairs of 
singly bound waters on opposite sides of the anion plane is 0.8 kcal/mol. 

The most stable SCF/DZP trihydrate complex Ilia has symmetry. 
However, a Djh optimization at the SCF/TZP+ level yields a double saddle 
point at -507.214732 hartrees. The degenerate imaginary vibrational frequen- 
cies lead to a minimum energy structure with a slight out-of-plane distortion 
of the complex (see Figure 3). The minimum energy structure lies only four 
microhartrees below the planar double saddle point structure. If pairwise 
ion-water and water-water interactions are additive, D e for the three-water 
complex Ilia is expected to be 3*14.3 - 3*1.2 = 39.3 kcal/mol. The complex 
is very slightly more bound than this ( D t — 39.4 kcal/mol). 


Atmospheric Implications 

Figure 8 gives the altitude profile of the mono- and dihydrate complexes 
as fractions of total nitrate considering only NO3 and H 2 0 in equilibrium, 
using a constant water concentration of 5 ppm by mass and the 1976 Stan- 
dard Atmosphere temperature and density profiles [35]. While quantitative 
predictions require consideration of other species (particularly HNO3), con- 
centrations of di- and trihydrate complexes in the stratosphere and lower 
mesosphere under these conditions are predicted to be very low. 
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Table 4: Theoretical versus experimental thermochemical quantities for the 
clustering of n water molecules on N0 3 . 


Method 


Source 



(Ha) 28.6 24.6 -25 

27.4 23.6 -24 

35.2 

(exper.) -28 


15.2 1 

14.3 

13.5 

18.5 
14.9 


13.1 -13.4 

12.3 -12.6 


-12.4 

-14.6 

-25.1 

-24.0 


(Ilia) | 40.5 34.9 -35 

(exper.) [ ^ -42 

e Boys counterpoise correcte 


-28.9 
34.9 -35.3 


-25.0 

-21.7 


-19.1 

-25.0 

-49.0 

-45.4 

-55.3 

-71.5 

-88.5 


SCF/DZP This work 

SCF/TZP+ This work 

SCF/QZ2P+ This work 

SCF/6-31G Ref. [12] 

SCF/DZP Ref. [13] 

HP MS Ref. [8] 

HP MS Ref. [6] 

SCF/DZP This work 

SCF/TZP+ This work 

SCF/6-31G Ref. [12] 

HP MS Ref. [6] 

SCF/DZP This work 

HP MS Ref. [6] 

or monohvdrate complex la are 


14.7 and 13.1 kcal for the DZP and TZP+ bases, respectively. 



Altitude (km) 
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Figure 8: Equilibrium distribution of nitrate hydrates in the stratosphere 
and lower ionosphere, assuming a water concentration of 5 ppmm. 



Equilibria under polar stratospheric conditions were modeled in the tem- 
perature range from 190°K to 200°K with an upper limit on the water con- 
centration equal to the vapor pressure of ice. 57% to 54% of the total nitrate 
present exists as monohydrate la and 36% to 40% exists as dihydrate Ila 
(considering only nitrate and water in equilibrium). The trihydrate Ilia and 
dihydrate lib each account for 2%. About 1% of total nitrate exists as the 
singly hydrogen bound monohydrate lb and only 0.6% of total nitrate is 
predicted to exist in unhydrated form under these conditions. 
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4 Conclusions 

Completely optimized equilibrium geometries, vibrational frequencies, ther- 
mochemical parameters, and stable site labilities of nitrate ion hydration 
complexes involving 1, 2, and 3 waters have been presented. The most sta- 
ble geometries in the first hydration shell involve in-plane waters bridging 
pairs of nitrate oxygens with two equal and bent hydrogen bonds. A second 
extremely labile local minimum involves out-of-plane waters with a single 
hydrogen bond and lies about 2 kcal/mol higher. 

The potential in the region of the second minimum is extremely flat and 
qualitatively sensitive to changes in the basis set; it does not correspond to a 
true equilibrium structure. A small but qualitative change in the geometry of 
the minimum energy trihydrate structure occurs in moving from the DZP to 
the TZP+ basis. Our results emphasize that the positions of mathematical 
minima on the energy hypersurfaces provide an incomplete description of the 
structure of weakly bound complexes. The potential around the minimum 
must be explored to characterize the lability of the complex, and this explo- 
ration should employ several different basis sets in order to detect basis set 
artifacts. 
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